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Abstract 

We compute the full 0(af) and ©(aglogOs) corrections to the potential VR(r) 
between the static color sources, where VR{r) is defined from the Wilson loop 
in a general representation i? of a general gauge group G. We find a violation 
of the Casimir scaling of the potential, for the first time, at 0{a'^). The effect 
of the Casimir scaling violation is predicted to reduce the tangent of VR{r)/CR 
proportionally to specific color factors dependent on R. We study the sizes of the 
Casimir scaling violation for various i?'s in the case G = SU{3). We find that 
they are well within the present bounds from lattice calculations, in the distance 
region where both perturbative and lattice computations of Vji{r) are valid. We 
also discuss how to test the Casimir scaling violating effect. 



PACS numbers: 12. 38. Aw, 12.38.Bx, 14.40.Pq 



1 Introduction 



The nature of the strong force is still a subject studied widely today. In particular, 
the static QCD potential has been studied extensively for the purpose of elucidating 
the nature of the interaction between static color sources. The static potential is a 
generalization of the Coulomb potential in QED to the case of QCD. Generally, the 
static potential at short-distances can be computed accurately by perturbative QCD. 
On the other hand, the static potential at long-distances should be determined by non- 
perturbative methods, such as lattice simulations or computations based on various 
models. 

Since some time, lattice computations have shown that the static QCD potentials 
between the color sources in various color representations (in the color-singlet channel) 
exhibit a property known as "Casimir scaling," within accuracies better than 5%, and in 
the distance range 0.1 fm < r < 1 fm lU |2j|j Casimir scaling is a property of the static 
potential Vr^t) between the color sources in the representation R, that the dependence of 
Vii{r) on R is given only by an overall factor Cr, the eigenvalue of the quadratic Casimir 
operator for the representation R. It turned out [5] that the Casimir scaling property 
of the static potential is a powerful discriminant of various models and approaches, 
which attempt to explain the nature of the QCD vacuum and color confinement. So far, 
however, a reasoning of Casimir scaling from the first principle has been missing. 

Computations of the static potentials in perturbative QCD have steadily made progress 
over decades. In particular, the discovery of the cancellation of the renormalons in the 
total energy of a static quark-antiquark pair led to a drastic improvement in the accuracy 
of the perturbative prediction of the static potential between the fundamental charges 
[B]. At the same time, the distance range, in which lattice computations and perturbative 
prediction agree with each other, extended to a significantly wider range [3 El [9] . By 
renormalization-group improvement of the perturbative prediction (after subtracting the 
renormalon), this overlap range extended to 0.05 fm < r < 0.4 fm [10]. It was shown 
that non-perturbative contributions to the static potential in this overlap region are 
much smaller than previous estimates, once the renormalon cancellation is incorporated 
in the perturbative computations. The linear potential (whose size is determined from 
the long-distance behavior of the potential) is excluded as a non-perturbative contribu- 
tion in this distance region. Furthermore, it was shown analytically that the perturbative 
static potential approaches a "Coulomb-|-linear" form in the above distance range due 
to the higher-order terms of the perturbative expansion, on the basis of a renormalon- 
dominance picture [TT] . 

Recently two groups (including our group) have independently completed computa- 
tions of the 3-loop corrections to the static potential between two fundamental color 
representations [12l[T3]. Compiling our present knowledge, the perturbative expansion 
of the static potential between the fundamental charges is known up to C(af) and also 
all the logarithmic terms at 0{al) [T3[[Tn] in this expansion are known. Comparatively, 
the perturbative expansion of the static potential between the color sources in a general 

^ See [2111] for simulation studies on the Casimir scaling property for gauge groups other than SU{3). 
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Figure 1: Contour C for the Wilson loop in 1^r[C]. 



color representation has been known only up to 0(a;^) [16j. The Casimir scaling is known 
to hold for the perturbative potential at least up to this order. 

In this paper we compute the full 0{a.'l) and loga^) corrections to the static 
potential between the color sources in a general color representation, for a general gauge 
group, by generalizing the computation for the potential for the fundamental represen- 
tation. The 0{a'^) correction originates from the purely perturbative 3- loop correction, 
while the 0{a'l\ogas) correction originates from the ultra-soft correction. We find that 
the Casimir scaling is violated, for the first time, at 0{a.^). In view of the consistency of 
the perturbative and lattice predictions for the fundamental potential, it is interesting 
to test consistency of the two predictions regarding the Casimir scaling. In particular, 
it is crucial whether the effects of the violation of the Casimir scaling by perturbative 
QCD exceed the current bounds from lattice computations. We examine the effects for 
various color representations within perturbative QCD and compare them with the lat- 
tice results. We also discuss how to test the perturbative prediction for the violation of 
the Casimir scaling. 

In Sec. 2, after fixing our notations and conventions, we present our result for the 
0{af) correction to the potential for a general color representation. The loga^) 
ultra-soft correction, as well as the complete expression for the potential up to 0{af) 
and C(a^logQ;s), are given in Sec. 3. Sec. 4 gives an analysis of the color factors in 
perturbative computations up to 3 loops. We examine the effects of the Casimir scaling 
violation for various representations in Sec. 5. Concluding remarks are given in Sec. 6. 

2 Wilson loop and 3-loop static potential 

We consider the (vector) gauge theory for a general gauge group G with n; massless 
fermions in the fundamental representation. In the case G = SU{3), the theory is QCD 
with rii flavors of massless quarks, which is our main concern. Nevertheless, to keep 
generality, all our formulas will be presented for the general gauge group. 

We define the static potential between the color sources in a general representation 
R, in the singlet channelo For this purpose we first define the vacuum expectation value 

^ Potentials in non-singlet channels are defined by inserting appropriate projection operators on the 
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Table 1: Some of the color factors related to the fundamental (vector) and adjoint represen- 
tations {R = F and R = A, respectively) in the case G = SU{N) and G = SO{N). Our 
convention is Tp = 1/2. 



of the Wilson loop as 



Wn[C] = (o|r 



Tr P exp l^zg, £ d^M;:(z)T^|] |0)/(0|Tr 1r\0) . (1) 



Here, T denotes the time-ordering of the gauge field operators ^^(-z). The hermitian 
generators of a general representation R are denoted by T^, which satisfy the commuta- 
tion relation [T^, T^] = if°-''^T^ with the structure constant f"'^'^ of the gauge group G. 
The fundamental and adjoint representations, respectively, will be denoted explicitly by 
R = F and R = A. The trace Tr is taken in the representation R, and Ir denotes the 
identity. P stands for the path-ordering of along the contour C, which is a rectangular 
loop of spatial extent r and time extent T; see Fig. [1] 
The normalization of is given by 

TY(T^T^) = TH^'^^ (2) 

We fix our convention by setting Tp = 1/2. The dimension of the representation R is 
denoted by Nr. The quadratic Casimir operator as well as symmetric invariant tensors 
which will appear in our computation are defined as 

= CnSij, (3) 

4-- = l^Tr(T^"«...T^"^"'), (4) 
n. ^ 

where the sum is over all permutations vr of the indices. There are a number of relations 
satisfied by these color factors; see [IB]. We note two of them: (i) TrNa = CrNr, which 



Wilson loop at temporal boundaries. Presently the potential for the fundamental representation in the 
octet channel is known up to 2-loop order fT7|. 
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follows from eq. ([3]) after taking trace; (ii) d'^'""'" = for odd n's, since (T^)'^ = —T%. 
Some values of the color factors are listed in Tab. [T]for G = SU{N) and G = SO{N). 

The static potential between the color sources in a general representation R is defined 
from the Wilson-loop expectation value as 

Vnir) = lim \og{Wn[C]) . (5) 

In other words, the potential appears in VTrJC] in an exponentiated form for large T, 
WR[C]r~^exp{-iTVR{r)}. 

We may evaluate Vji{r) using perturbation theory. Since our technology for loop 
computations is developed mostly in momentum space, we compute the potential in 
momentum space. In order to regularize both UV and IR divergences, we employ di- 
mensional regularization with one temporal dimension and d = D — 1 = 3 — 2e spatial 
dimensions. Thus, the perturbative potential expressed in terms of the corresponding 
(l^-scheme) coupling in momentum space reads 



(6) 



where both quantities are denoted with superscripts PT to make explicit that they are 
computed in perturbative expansions of the strong coupling constant. A prefactor is 
included on the right-hand side such that ayj(g) is defined to be dimensionless; g = 
je = 0.5772... denotes the Euler constant. 

The perturbative expansion of ayj((?) is expressed as 



oo 



a 



■7j 



n=0 

with 




(8) 

Here, as(/i) denotes the strong coupling constant renormalized at the renormalization 
scale /i, defined in the modified minimal subtraction (MS) scheme; Pn{L) denotes an 
rath-degree polynomial of L. The renormalization-group equation of as{^i) is given by 

rflog(/.2) 1, An ) [ An ) ' 

where /3„ represents the {n + l)-loop coefficient of the beta function. The relevant 
coefficients read [19] 

11 ^ 4 _ ^ 34 ^2 /20 



/3-i = e, /3o = yC^A-3ri,T^, /3i = — Ci - Ca + ACp j riiTp , 

158 ^ 44 \ 2 2 nn^ 
— CA + —CFjniTp. (10) 
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Both V^^{r) and ay^{q) are independent of the scale /x, hence their dependences on 
log yU^ are dictated by the renormalization-group equation. For n < 2, the only part of 
the polynomial Pn{L) that is not determined by the renormalization-group equation is 
an = -Rrt(O). For n > 3, Pn{L) includes IR divergences in terms of poles of e and assoicated 
logarithms, whose coefficients are not determined by /3j's. Up to 3-loop order, they are 
given by 

Po = ao. Pi = ai + aoPoL, P2 = 02 + (2ai/3o + aol3i)L + ao/So^L'^, 



Ps = as + (3a2/3o + 2ai/3i + ao/Js)^^ + ( 3ai/3o' + -aoMi ]L' + aoPo L 



3 rS 



and 



as 



+ 3L 



(11) 



(12) 



The 1/e term represents the IR divergence [20 
the corresponding scale dependence. 

The first three independent of the representation R 



, which is separated together with 
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(13) 



where Cs = C(3) = 1.2020... denotes the Riemann zeta function C,{z) = J2'^^=i^/'n^ 
evaluated at z = 3. Since Cr is factored out in eq. (jH]), and since /3„'s are independent of 
the representation R, if all a„'s are also independent of R, the potential V^^(r) satisfies 
the Casimir scaling. 

The 3-loop non-logarithmic constant 03 for the fundamental color sources has been 
computed recently: the contributions with internal fermion loops in [23] and the purely 
gluonic contributions in [T^l US]- We repeat these calculations for the general color 
sources. We find 



as = - 
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■ / 12541 647r^ 
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(14) 



where Cs = C(5) = 1.0369.... In the above equation, the coefficients Cj's are known only 
numerically: 



Cl = -354.859, C2 = -56.83(1), 



(15) 
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C3 = 502.22(12) , C4 = -136.8(14) . (16) 

We find that 03 is dependent on the representation R through the color factors d'^'^'^d'^'^'^/ {NaTh) 
and d°^'^'^d'^^'^'^ / {NaTji). These are the leading effects of the Casimir scaling violation in 
perturbation theory. 

One may want to include massless fermions in representations other than the fun- 
damental representation. One may easily identify the contributions from fermions in a 
general representation R' by simple replacements Tp — )■ Tr/, Cp Cri and d'^p'^'^ — d'^^F'^ 
in the coefficients of rii in the equations presented in this section. 



3 Ultra-soft correction 

The IR divergence contained in 03 is an artifact of the strict perturbative expansion of 
the static potential VR(r) in as- The divergence originates from the degeneracy of the 
singlet and adjoint intermediate states in PVii;[C] within the naive perturbation theory 
|20j . Beyond (naive) perturbation theory, this degeneracy is known to be lifted and 
the IR divergence is absent. The difference between the static potential VR(r) and its 
perturbative expansion V][^{r) can be treated systematically within the effective field 
theory Potential Non-relativistic QCD (pNRQCD)0 [21], which treats ultra-soft (US) 
gluons and heavy quark-antiquark bound-states as dynamical degrees of freedom. This 
difference VR(r) — V^'^(r) is given by contributions of ultra-soft degrees of freedom. 

We repeat the computation of the ultra-soft contributions given in [211 125] for the case 
of VR(r), taking into account the following points: (1) We change the representation of 
the static charge of the Wilson loop from the fundamental to the general representation 
R. (2) The perturbative matching between the full theory and pNRQCD is properly 
taken into account, such that the non-logarithmic term of VR{r) — V^^{r) is predicted 
correctly at 0{af). Thus, we have 



5V^\r)^Vn{r)-Vr{r) 



Jo 



-Vs) 



x(o|r 



r-E"(t)fPe*^^/o^*'^S(OTl\ r-E''(0) 

V J ab 



\0) + O{r^). (17) 



The term shown explicitly is the leading-order [(9(r^)] term in the multipole expansion 
in r, where r = |r|. All fields in the matrix element are evaluated at the spatial origin. 
{T^)ab = —if"''"^ denotes the generator of the adjoint representation. Vg and Kd denote 



3 The results of Ref. [13] are C3 = 502.24(1), C4 = -136.39(12). 
Strictly speaking, the name of the effective theory, as well as terminology such as "quarks" and 
"gluons", are adequate only in the case G — SU{3). We will frequently use them for the general gauge 
group G, nonetheless. 
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Figure 2: Class of diagrams contributing to [<JV^^]lo- Dashed lines represent Coulomb gluons, 
while a curly line represents transverse gluon, in Coulomb gauge. 



the singlet and adjoint potentials, respectively, in the representation i?, i.e., the potentials 
between the sources in the singlet and adjoint channels, respectively, in i? ® -R. 

In the distance region r -C Aq^j^, we may further expand the matrix element as well 
as Vad — K in the coupling constant and obtain the leading-order contribution to SV^^^r) 
in double expansion in and log(Q;s)|2ll§ It is understood that we employ dimensional 
regularization, and Qs can be identified with the bare coupling constant of the full theory 
in the leading order. Noting that the singlet and adjoint potentials are given by 



VJr) 



Kd(r) 



ifi'^r^e^^Y r(l-2e) 



rfi 



+ 0{a 



Sj 5 
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S) 5 



(18) 
(19) 



we obtain 



LO 



-gl Cr 



d-l 

d 



■(K 



ad 



47r 



T{l + d)T{-d) 



2A-nr 



(20) 



where we dropped (9(e) terms in the last expression. Since [5V^^]lo depends on the 
representation R only through the overall factor Cr, Casimir scaling is preserved by this 
contribution. The 1/e term represents the UV divergence of [5V^^]lo within pNRQCD 
effective theory. Upon Fourier transform, the 1/e and log(/i^) terms of eqs. (fT2|) and (120|) 
cancel each other. 

We made a cross check of eq. (120]) as follows. It was first pointed out in [20] that 



[(JV^^Ilo may be obtained as the difference between the resummation of Coulomb ladder 



^ One can verify that indeed 5V^^{r) = ^^(r) — y^"'"(r) vanishes in the strict perturbative expansion, 
by expanding also e~** in a^; in this case, each integral becomes scaleless, and therefore it 
vanishes in dimensional regularization. 

^ 0{r^) terms in eq. ([T7|) do not contribute to the leading order of this double expansion. 
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diagrams, given in Fig. |2l and the sum of the individual diagrams (in Coulomb gauge) 13 
Namely, [5V^^]lo inay be computed more directly using a resummation of diagrams in 
perturbation theory, without recourse to pNRQCD effective theory. One can verify that 
this is indeed the case, including the non- logarithmic terms in eq. f l2U]) . Formally the 
equivalence of the two approaches can be shown using the threshold expansion of Feyn- 
man diagrams [26] . We have also evaluated the relevant diagrams directly and confirmed 
the coefficients of the logarithms as well as the non-logarithmic constant explicitly. In 
the resummation approach, one can avoid complexity related to the matching of the full 
theory to pNRQCD, in computing the non- logarithmic terms. (See also |27j.) 

For completeness, we present the formula for VR(r), as given by the sum of V^'^lr) 
and 6V^^{r), including all the corrections up to O^a^) and 0(a;^ logo;,,). Namely we use 
the series expansion ([7]) up to = 3 for the former and eq. fl20|) for the latter. We obtain 

>'«wL™.o-^«^i:/"(^')(^)". (21) 

where 

L' = log(/iV) + 27£;, (22) 

and 

pQ = ao. Pi = ai + ao(3oL', = a2 + (2ai/3o + ao(3i)L' + ao/3o^ (l''^ + ^) , 
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Pg = as + 5a^s + (Sag/^o + 2ai/3i + aof32)L' 
5 



+ (3ai/3o' + -ao/3o/3i) + y) + ^oPo'{L" + vr^L' + I6C3) , (23) 



with 



3 



log(C^as(/i)) +-fE- - 



(24) 



By definition [VR(r)]NNNLO is a renormalization-group invariant quantity and is free from 
IR divergences. Apart from the overall factor Cr, only in the above expression is 
dependent on the representation R, which is the source of the violation of the Casimir 
scaling. 

Let us comment on the non-perturbative contributions to the static potential Vii{r). 
At r < Aq^j-,, the multipole expansion in r of 5V^^{r), given in eq. f|T7|) . constitutes an 
operator-product-expansion of 6V^^{r). In this way, one can systematically parametrize 
the non-perturbative contributions to SV^^{r) [and therefore to VR(r)] in terms of ma- 
trix elements of non-local gluon condensates, after subtracting UV divergences. See 
PT| O [ini ES] for analyses of the non-perturbative contributions to V/j(r) in this frame- 
work. [5V^^]lo can be regarded as defining a scheme for subtracting the UV diver- 
gence at leading order, or equivalently, as defining a scheme for the corresponding non- 
perturbative contributions. 



^ To be precise, one should replace Ca/2 — Cr by Cyi/2 in the latter (individual) diagrams, in order 
to take into account the shift of the ground-state energy ^> 14 (r) by Coulomb resummation. 
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Figure 3: Class of diagrams contributing to the potential at 1-loop in Coulomb gauge. The 
hatched part represents the vacuum polarization of the Coulomb gluon. 

4 Color factors in loop diagrams 

We examine the general structure, especially the R dependence, of the color factors in 
V^^{r) without going into details of the loop integrals. A discussion on this subject up 
to 2 loops is given briefly in [16J. We extend the argument to 3 loops and to more details. 

At tree-level, the color factor of the potential is simply Tr(T|^T^)/Tr(l/j) = Cji. 
As is well-known, a simple way to understand the color factor at 1-loop order is to 
consider the corrections in Coulomb gauge. In this gauge, only the vacuum polarization 
of the Coulomb gluon contributes to the potential [29]; see Fig. [31 Since the vacuum 
polarization is independent of R and is proportional to 5"^, the i?-dependence of the 
potential is the same as the tree-level potential. The same argument cannot be used for 
the 2-loop or higher-order corrections, however, since there are contributions from graphs 
other than the vacuum-polarization type. Hence, we develop a more general argument. 

Let us make some preparations. First, according to eq. Tr(T^^ ■ ■ -T^") equals 
'^fl plus a sum of the terms which include at least one commutator [T^%T^^]. If 
we rewrite the commutator as if°-i°-i^T^ and repeat this operation recursively, we may 
express the trace by d°^"'""^ and lower rank tensors as 

n-1 

Tr(T^^ ■ ■ ■ T^^-) = 4-"" + dn{k) , (25) 

k=2 

where in the second term on the right-hand side we suppressed the color indices; dR{k) 
represents a fcth-rank invariant tensor; c is an i?-independent coefficient and expressed 
in terms of f"'^'^. {dfi{k) and c have appropriate color indices which are not shown.) 
Contracting both sides with f°-i°-3^^ where {oj, aj} C {ai, . . . , a„}, the first term on the 
right-hand side drops. Hence, we have 

n-1 

pa,b T^^yai . . . T^") = ^ ^ c' ■ dR{k) , (26) 

k=2 

where we absorbed in c' . Note that on the left-hand side there are n generators 

T^, while the highest rank of the invariant tensor on the right-hand side is — 1. 

In perturbative computations of the static potential, one first eliminates iterations 
of lower-order potentials from the diagrams contributing to the expectation value of the 
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Figure 4: Reduction of a 2PI color graph to a 2PR graph and a residual graph using the 
identity T^T^ = Tj^T^ + if^^T^. The 2PR graph (first graph on the right-hand side) will be 
deleted. In the second graph, the number of vertices on the static line is reduced by one. 



Wilson loop, VI/ij[C]. A prescription has been developed, in association with a proof 
for the exponentiation of the static potential [30] • Let us briefly review the prescription. 
Within the diagrams for iyR[C], only the 2-particle-irreducible (2PI) diagrams contribute 
to the potential. Here, 2-particle irreducibility is defined with respect to cutting the two 
static time-like lines simultaneously. Furthermore, among the 2PI diagrams, the color 
factors of the diagrams, which become 2-particle reducible (2PR) by sliding vertices on 
the two static lines, need to be modified; this is because, part of these diagrams are 
identified with iterations of lower-order potential]^ How to modify the color factors 
of these diagrams is the essence of the prescription. One associates a "color graph" 
with each of these diagrams. Let us denote the original 2PI diagram by D and the 
corresponding color graph by Cq{D). To begin with, the topology of the color graph 
Cq{D) is taken to be the same as that of D, and the value of Co(-D) is set equal to the 
color factor of the diagram D. Then one tries to render the color graph Co{D) to a 2PR 
graph by sliding vertices on the two static lines. Whenever any two vertices a and 6 on a 
same line are interchanged in this procedure, one separates Co{D) into two color graphs 
Ci{D) and 02(0), corresponding to TgT^ = T|T^ + [Tg,T^]. In the latter graph, the 
two vertices are combined to a single vertex on the line, corresponding to i/^'^^T^; see 
Fig. m Thus, the value of C2{D) equals the value of Co{D) except that T^T^ is replaced 
by if°''"^T^. This procedure is repeated for each color graph, until no more reduction of 
any of the 2PI graphs into a 2PR graph is possible. Then one deletes all the 2PR graphs. 
The sum of the values of the remaining 2PI color graphs I]je2Pi Ci{D) is the color factor 
that should be assigned to the original diagram D. 

After applying the above prescription, (the value of) any of the color graphs Ci{D), 
which correspond to a general 2PI diagram D, can be written in the form 



= NrCr X 



(27) 



Na NaTr 

where X^^ ' ^^'s are independent of R. In the second equality, we substituted eq. (|25 



^ If the gauge group is abelian, any diagram, which becomes 2PR by sliding vertices on the two 
time-like lines, does not contribute to the potential. This is because, these diagrams (at a fixed order) 
sum up exactly to the exponentiation of the lower-order diagrams which contribute to the potential, if 
we can ignore orderings of T^'s. 



10 



(b) 

Figure 5: (a) 2-loop diagram D with n = 6, and (b) the corresponding color graph C{D) with 
Uc = 4. Starting from the diagram (a), move the right-most vertex on the upper Hne to the 
left by interchanging the vertices twice; after removing the 2PR graphs, one is left with the 
2PI color graph (b). 



a 



and used d!^^ = Tr5'^*, Tr = NfiCr/Na- The overall factor Nr will be canceled in the 
potential by the denominator Tr(lii;) = Nr. Each term of x'^^'""''' is not factorizable (e.g. 
j^abcd ^Qgg j^Q^ include a term of the form A°'^B'^'^), since a factorizable term corresponds 
to a 2PR color graph@ We will see this in explicit examples below. 

We note that, in eq. fl27|) . purely gluonic contributions to - (expressed in terms 
of only /"^'^'s) for an odd k vanish, at least up to fairly high orders in perturbative 
computations, and it is definitely true up to 3 loops. This is because, X'^^'""''' can be 
expressed with only invariant tensors c?^^'"^^'s (/'^'"^'s can be eliminated) up to certain 
high orders [18]. Since cf^'""'^ = for an odd j, j^ai---ak gj^Q^y vanish!^ 

Now we are in a position to discuss the loop diagrams. We start from 2-loops. Let 
us denote by n the number of vertices on the Wilson loop in a 2PI diagram D. Then 
we examine, for each fixed n(< 6), the corresponding color graphs. Consider a diagram 
D with n = 6 or n = 5. The diagram does not include fermion loops. The number of 
vertices on the Wilson loop in the corresponding color graphs C (D) is tIc = 4; see Fig. [51 
The color graphs have a form C{D) = Tr(T^i ■ ■ ■ T^*) ■ Each term of X^^-'^* is 
expressed by two structure constants J^^^ jcde ^j^^j-^ different assignment of indices. We 
rewrite Tt{T^^ ■ ■ ■T'^) in terms of cLrs using eq. (125|) . Then C{D) is expressed in the 
form of eq. (127]) . There is, however, no way of contracting indices such that d'^^'^''' survives. 
Since Xf^'^ is zero in the pure gluonic case, the only remaining term is Xf". 

We consider the diagrams D with n = 4. The corresponding color graphs have 
either of the following 3 forms: (i) C{D) = Tr(T^i ■ ■ ■ T^*) with each term of 

Xf expressed by two structure constants; (ii) C{D) = Tr(T^T^T^) Xg"^^; (iii) C{D) = 
Tr(T^T^) Xf''. The case (i) is the same as n = 5,6. In the case (ii) the purely gluonic 

^ 2PR color graphs, which are deleted, contain powers of Cu if they contain (for instance) iterations 
of the tree graph, X'^^--'"" oc J'^i"^ ^"3^4 . . gjnce T^Tg = Cr1r._ 

10 Furthermore, if we take a sum Q{D) + C,{D), where C-{D) = (-l)"Tr(rj^" • • • Tj^^) X^i- '^", in 
the sum, purely gluonic contributions to ^'^ - ''-k x^^- -^'' for an odd k vanish to all orders. Note that, 
due to the charge conjugation symmetry, if there is a color graph Ci{D), there is also a conjugate graph 
Ci{D), unless the two graphs are the same. 
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(a) (b) (c) (d) 



Figure 6: Pairs of color graphs with a fermion loop: (a)(b) 2-loop graphs; (c)(d) 3-loop graphs. 
The graphs in each pair have the same topology and the opposite charge flow directions in the 
fermion loops. The d'p'^ part of Tv{TpTpTp) cancel between each pair. 

contribution to the d'^'^Xf"^ term vanishes, as we already noted; on the other hand, the 
contributions of graphs with a fermion loop to d'^'^Xf"^ is proportional to /"'"^d^^'^ = 0, 
since the fermion loop is included as a gluon vacuum polarization. The case (iii) is of 
the type by itself. Thus, C{D) has only the term in all three cases. 

The color graphs, which corresponds to the diagrams with n = 3, are either in the 
form Tr(T^TjjT^) X!^'"^ or Tt{T^T^) Xf'. The only new aspect, as compared to the cases 
(ii) and (iii) of n = 4, is in the diagrams including a fermion loop which has three 
vertices on it. These diagrams are proportional to n,Tr(r^T|T^) = rii^df" + fT^/"^^). 
The nid'^'^ part is canceled between a pair of color graphs with the same topology but 
with the opposite charge flows of the fermion loop; see Figs. [6]|^a)(b). This is because 
TT{TpTpTp) — TiiTpTpTp) = iTpf"'^'^. Thus, the sum of such a pair of color graphs is 
niTp times the corresponding 1-loop color graph, in which the fermion loop is replaced 
by the gluon three-point vertex. 

The contributions of the diagrams with n = 2 is of the type X|". Hence, all the 
contributions are reduced to Xf" of eq. ( 127|) . which is independent of R apart from the 
overall factor N^Cr. Thus, the Casimir scaling holds. 

We can discuss the 3-loop case in a similar manner. As before, we discuss the color 
factors of 2PI diagrams for each fixed n, the number of vertices on the Wilson loop. The 
maximum number is n = 8. 

For a moment, we consider diagrams without internal fermion loops. In the cases 
n = 8, 7, 6 and 5, the number of vertices on the Wilson loop in the corresponding 
color graphs C{D) is ric = 5 or less. In the case ric = 5, the color graphs have a 
form C{D) = Tt{T^' ■ ■ ■ T^^) X^'-"^, and each term of X^'-"^ contains f^^^'' ({a^, aj} C 
{oi, . . . , fln})- According to eq. ( l26l) . the highest rank invariant tensor in C{D) is the 4th- 
rank tensor d'^^'^'^. C{D) includes the Xg" and X|*'^'^ terms of eq. ( 127|) . Each term of X^'^'^^ 
includes four f^^'^ = i(T^)gh and can be expressed a£3 Tr (T^'T^'T^'T^') with a different 
assignment of (a', b', c', d') to (a, b, c, d). There are no contractions of the indices of f^d^^s 
in the form Ti[T% T\ )Ti{T% ) oc 5" ; these correspond to 2PR color graphs and 
would have generated C\ had they contributed. Using eq. ( I25l) for R = A, d'^^'^'^Xf"^'^ can 
be expressed by d^^^'^df^'^. Thus, C{D) is expressed by X^'^ oc C^Na and Xf^'^ oc df""^. 

The arguments in the other cases {n = 2, 3, 4 and ric < 4 for n = 5, 6, 7, 8) are the 

Many contractions are trivially zero, such as f^be^abcd _ -^^Jijch are not discussed explicitly. 
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Figure 7: Some of the 2PI diagrams which contribute to d'^^'^'^d'p^ and df^^'^'^d"^"^ . 

same as above, apart from the point that for n = 2, 3 there is no Xf"^'^ term. In the end 
the color factors are expressed by oc C\Na and Xf'^'^ oc d"^'"''. 

Now we turn to the diagrams with fermion loops. We denote the number of vertices on 
a fermion loop by n'(< 6). A fermion loop with n' = 2 is proportional to niTi{TpTp) = 
riiTpS""^. Hence, if a 3-loop color graph includes at least one fermion loop with n' = 2, 
it is equal to riiTp times the corresponding 2-loop graph, in which the fermion loop 
insertion is removed. In this case, the R dependence is only in the overall factor NjiCr. 

A fermion loop with n' = 3 is proportional to niIi{TpTpTp) . As in the 2-loop case, 
the nid'^p'^ part is canceled between a pair of color graphs with the same topology but 
with the opposite charge flows of the fermion loop; see Figs. IS^c)(d). Thus, the sum of 
such a pair of color graphs is n/Tp times the corresponding 2-loop color graph, in which 
the fermion loop is replaced by the gluon three-point vertex. 

A fermion loop with n' = 4 is proportional to niTY{TpTpTpTp) , ni'Ti{TpTpTpTp) or 
uiTiiTpTpTpTp) , and a 3-loop color graph includes one such loop at most. The latter 
two traces are proportional to riiTpd""^ and represent gluon vacuum polarizations; namely 
the color factors in these cases are proportional to those of lower-loop graphs. The first 
trace generates nid'^'^'^d'^'^'^ and lower rank tensors via eq. (l25l) . Also in this case, the 
terms with nid'^'^ are canceled between a pair of color graphs with the same topology 
but with the opposite charge flows of the fermion loop. This can be seen as follows. Due 
to the hermiticity of T^, Tr(T^T|T|T|5) + Tr(T^T|T^T^) = 2 Re [Tr (T^T^T|T|5)]. If we 
substitute eq. f l2^ . the terms with odd numbers of f^^^ are pure imaginary and vanish, 
hence only the even rank tensors d'^'^'^ and dp = TpS""^ remain. 

A fermion loop with n' = 5 is proportional to Tr(Tp^ ■ ■ ■ T'p) and included in diagrams 
with n = 2 and 3. In most of the corresponding color graphs, at least one pair from the 
indices (ai, . . . , 05) is contracted between the pair; in this case, the trace can be reduced 
to TiiTpTpTp), and the riid'^p^ part is canceled between a pair of color graphs with the 
opposite charge flows of the fermion loop, as before. The only other type of contraction 
is in the form Tr(T^T^T|T^T|)d^''/"^^ which can be reduced to d^", i.e. the X!^'' term. 

A fermion loop with n' = 6 is contained only as a gluon vacuum polarization in the 
diagrams with n = 2. Hence, they contribute only to the term of eq. f l27|) . 

In summary, up to 3-loops the Casimir scaling violating terms arise only as d'^'^'^d'^'^'^ 
and nidf^'^d"^"'^ in the Xf"'^ term of eq. (EZD- Typical 2PI dia grams which contribute 
to these color factors are shown in Fig. [71 All the other color factors are collected in 
the term. The absence of the odd-rank symmetric invariant tensors dpidA follows 
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from (T^)"^ = —T^ and the cancellation between diagrams with opposite charge flows in 
a fermion loop. These features are ensured by the charge conjugation symmetry. Thus, 
the charge conjugation symmetry plays an important role in suppressing the Casimir 
scaling violating effects. 



5 Numerical analysis of Casimir scaling violation 

In this section, we examine the perturbative predictions for the Casimir scaling violation 
numerically and compare them with lattice computations, in the case of QCD with zero 
flavor of quarks (quenched approximation). For the lattice results, we will use those 
given in [2], as they seem to be most accurate and extensive up to date. Throughout 
the analysis, we use the central value of roA^°^ = 0.574 ± 0.042 |I0] to fix the relation 

between the lattice scale and A|^°'', where tq denotes the Sommer scale j3l]. (It is 
customary to interpret tq ~ 0.5 fm when comparing this scale to one of the real world.) 
When making comparisons, one should note the following point. As is well known, the 
constant (r- independent) part of the potentials by perturbative and lattice computations 
cannot be related unambiguously. This is because, the self-energy contribution and 
the potential energy contribution in a lattice computation of the Wilson loop cannot 
be separated unambiguously. Furthermore, it is formidable in perturbation theory to 
reproduce exactly the subtraction scheme used in [2], which subtracts the leading-order 
self-energy contribution defined in lattice perturbation theory. Alternatively, we study 
quantities which are, by construction, independent of the constant part of the potentials, 
and discuss their relation to the Casimir-scaling violating effects measured in the lattice 
simulation. 

In this analysis we define the following two quantities as measures for violation of the 
Casimir scaling of the static potential for the general representation R: 

_ [VR{r)-VR{n)]/Cn ^ 

- [V,{r) - V,{n)]/C, ^'^^ 

Both quantities are zero if the potential has an exact Casimir scaling property. Both 
quantities are independent of the constant part of the potential: In the first quantity 
Sqq, we subtracted a constant from the potential VR(r) such that it vanishes at r = ri; 
in the second quantity 5^g, we used the force Ffi{r) = — V^(r) between the color sources. 
The fundamental representation {R = F) is used as a reference in both quantities. 

In principle and can be measured in lattice simulations. On the other hand, 
they can be computed in perturbation theory, by replacing VR,ir(r) by [Vr,f('")]nnnlo- 
They are equal at the leading order: 



C 71 / ^CL^cd^f^bcd JabcdJabcd'' 



F "-R 



2Na \ Tr Tf 
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+ 



C4 

2Na 



' ^abcd^abcd 



R 



^abcd^abcd ^ 

TV 



0{a 



(30) 



At this order, and are independent of r; is also independent of ri. 

Eq. (EnD shows that and 5cf are negative and their absolute magnitudes increase 



with the color factor rii d 



''^bcd^abcd^l^pj^j.^-^ or with df^df'''^/{NATR). (Note that C2, C4 < 
0.) As a result, the effect of Casimir scaling violation is to reduce the tangent of Vuir) / Cr 
for larger rii d'^p^'^d"^'^ / {NaTr) or for larger d'^^'^d'^^^/ (NaTr). This can be seen from the 
fact that Sqq, defined from the static force, is reduced by this effect. 

As shown in [12], in the case G = SU{3), R = F and rii = 0, the potential 
[Vp(r)]NNNLO for the fundamental representation agrees fairly well with the lattice results 
in the distance range 0.1 < r/ro < 0.50 Hence, it would be natural to test whether the 
perturbative prediction is also capable of reproducing the Casimir-scaling violating ef- 
fects in the same distance range. Generally the perturbative predictions for and 6q^, 
eq. (1301) . are quite small, since they are 3-loop effects. At the same time, they turn out 
to be strongly scale-dependent, since only the leading-order terms with a high power of 
as are known. For example, for the adjoint (octet) potential of SU{3) {R = A), we find 

■^^ -0.129 - 0.0030 rii) 



,Va ^Fa 



0.00032 (ni = 0, /i = /ii = Ag°P/0.035 



aS°7o.i4 



(31) 



-0.0013 [ni =0, /i = /i2 

Here, the scale Hi {fj.2) is the largest (smallest) scale0 with which the stability of 
[^('^)]nnnlo has been examined in [12j, where /^i and /i2 differ by a factor of 4. It 
is a standard way to estimate uncertainties of the perturbative prediction for the poten- 
tial. Since and are determined by the potential, we use this scale dependence 
to estimate uncertainties of these quantities. Hence, our estimates for ^Qg and ^^s 
between —0.03% and —0.13% for ni = 0, in the distance range 0.1 < r/ro ~ 0.5. The 
large uncertainties in the estimates are due to the strong scale dependence. The scale 
dependence is expected to be reduced if the higher-order corrections (beyond 3-loop) to 
5^3, (5cs incorporated. 

We list the color factors for various representations of SU{3) in Tab. [21 Our estimates 
for ^Qg, for the corresponding representations are shown in the same table for rii = 0. 
All the estimates are derived in the same way as in eq. fIST]) . Note that in the case of 
SU (3) simple relations 



JabcdMibcd ' 
F R 






NaTr 


SU{3) 




^abcd^abcd ' 






NaTr 


SU{3) 





Cr 



Cr 



SU{3) 



5(7(3) 



(32) 



(33) 



It is estimated in Ref. [TU] that, by renormalization-group improvement of [VR(r)]NNNLO (after 
subtraction of the renormalon), the agreement becomes better up to larger distances, r/ro ^ 0.8. 



They correspond to /^^ ^ « 0.06 ro and 1^2^ ^ 0.24 tq. 
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Nr 


( T) (7 1 






^abcd^abcd 
F R 


^abcd ^abcd 


5^5 6!^§(%) for 727 = 

CS' CSV'^/ ' 




NaTb. 


NaTr, 


3 


fl 0^ 


4/3 


1 /2 


5/48 


1 5 /8 




8 




O 




5/16 


45/8 


_n (13 1 3 


6 


f2 0) 


10/3 


5/2 


17/48 


51/8 


_0 04 16 


15a 


(2,1) 


16/3 


10 


29/48 


87/8 


-0.08 - -0.31 


10 


(3,0) 


6 


15/2 


11/16 


99/8 


-0.09 - -0.37 


27 


(2,2) 


8 


27 


15/16 


135/8 


-0.13 - -0.52 


24 


(3,1) 


25/3 


25 


47/48 


141/8 


-0.13 - -0.55 


15. 


(4,0) 


28/3 


35/2 


53/48 


159/8 


-0.15 - -0.63 




■(fc+l)(fc+2)' 
2 


{k,0) 

s 


fc(fc+3) 
3 


fe(A:+l)(fc+2)(fc+3) 
48 


2fc2+6fc~3 
48 


3(2fc2+6fc-3) 

8 


(-0.0063 - -0.026) 
x(fc+4){fc-l) 



Table 2: Color 
3, 8, 6, 15a, 10, 27, 24, 15,, 



factors 

(fc+l)(fc+2) 



for SU{3), 
, where each 



for the representations R = 
R is labeled by its dimension Nr. 
{p, q) denotes the weight factor. See Sec. 2 for the definitions of the other color factors. 
The last column represents the estimates (in per cent) of '^CS''^CS ^'^'^ ~ ^- '^^^ 
represents the A:th-rank completely symmetric representation. 



hold, since [df''%u{^) = [df^'^]su(3)/^8 = {5^''S'''^ + S'^'^S'"^ + S'^'^S''^) /2A. It follows that the 
contribution of the internal quark per each flavor is about 40 times smaller than that 
of the pure gluons, independently of R. Hence, the contributions of the internal quarks 
may be neglected in QCD in the first approximation. Eqs. ( |32l) and (l33l) also show that 
the Casimir-scaling violations, and 6q§, scale proportionally to Cr for Cr ^ 1 in 
QCD. 

Let us compare our estimates with the lattice results. We may relate to the 
Casimir scaling violating effects measured in the lattice simulation. From Fig. 4 and 
the numbers hsted in Tables II and IX of [2], Roir) defined in eq. (31) of that paper is 
estimated to be close to {Cr/Cf){)- + ^cs) with a choice ri ~ 0.05 tq. Or equivalently, 

^ll ^ ^Rd - 1. (34) 

Since our prediction for Sq^ is independent of ri at leading order, we can compare our 
predictions with the lattice results without a precise knowledge of ri. Ref. [2] observed 
no violations of the Casimir scaling in Rd at 0.3 ^ r/ro ^ 2 within 5% accuracy, for 
the representations R = 8, 6, 15a, 10, 27, 24, 15,. This translates to bounds < 5% 
for all of these representations in the same distance range. As seen in Tab. El these 
bounds are perfectly consistent with our predictions in the overlapping distance range 
0.3 ^ r/ro ^ 0.5: our estimates are an order of magnitude smaller in the largest cases. 

It is desirable to test our predictions in lattice simulations. On the one hand, our 
current estimates for Sq^, have large uncertainties due to scale dependences, and 
only their orders of magnitude can be predicted. On the other hand, the prediction, that 
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these quantities scale proportionally to specific combinations of color factors as given in 
eq. fl30|) . is expected to be more secure, especially for smaller r. As already stated, the 
Casimir scaling violation is enhanced for larger Cpt in QCD. Hence, if a representation 
with an extremely large dimension is chosen, the violation may be detectable. For 
instance, if we use the fcth-rank completely symmetric representation with k = 30, the 
estimates for ^^g, amount to between —6% and —26%; see Tab. [2l 

Another possible direction is to consider larger groups. As can be seen in Tab. [H 
the color factor df^'^df^'^ / {NaTa) increases rapidly with N both for G = SU{N) and G = 
SO{N). However, the terms with the leading power in N cancel between d'A'^'^d'A'^'^ / {NaT a) 
and df^'^dp'^/{NATF). As a result 5^^, S^i scale only linearly on N for SU{N) and 
quadratically for SO{N). Thus, even if a relatively large N is chosen, it seems necessary 
to use higher rank representations in order to detect Casimir scaling violating effect. 

There is a way to perform an indirect test, presumably with less efforts. Although 
the leading N dependence for G = SU{N) or G = SO{N) is canceled in 6q^, 6q^, this 
dependence is included in the fundamental potential Vp^r). In principle, it is possible 
to test the contribution of the d^^'^'^d'^p'^'^ / {N aTf) term in Vp{r) by comparing lattice 
computations and [Vf('")]nnnlo- In the case G = SU{N), the contribution proportional 
to d^^'^'^d'^p"^ / {N aTf) in [VF(r)]NNNLO (or in [dVpir) / dr]^^^i^o) is estimated to be between 
—6% and —23% of the total amount for = 25, and between —10% and —39% for 

= 30, using the same values for and as(/i2) as the SU (3) case. We may also use 

larger representations and test the d^^'^'^d'^^'^'^ /{NaTr} term. Verifying the contribution 
of the d°A'^'^d'^^'^'^ / [NaTr) term may be a first step toward testing the Casimir scaling 
violation. 

6 Concluding remarks 

In this paper, we computed, within perturbation theory, the full 0{af) and (^(a^logas) 
corrections to the potential VR{r) between the static color sources, defined from the 
Wilson loop in a general representation i? of a general gauge group G. The strictly 
perturbative contributions and the ultra-soft contributions are separately presented. The 
full expression for the potential up to 0{a^) and loga^) is also given. 

One way to utilize our result is as follows. In the case G = SU{3), the present result 
provides many observables (corresponding to different i?'s) which can be used in the 
matching procedure between lattice QCD and perturbative QCD. Already there exist a 
number of computations of the potentials for various i?'s by lattice simulations. Since the 
matching between both theories using the fundamental potential has been successful and 
can be used for the determination of as [IQ], it is expected that our present result may 
be used in the same way and may contribute to improve accuracies in the determination 
of or in predicting other physical observables. 

A prominent feature of our present result is that the violation of the Casimir scaling 
of the potential is predicted at (9(a^). This is the first prediction, based on a first- 
principle calculation, of the Casimir scaling violation for the QCD potential. The effect 
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of the Casimir scaling violation is to reduce the tangent of Vpi{r)/CR proportionally 
to df^'^df'-'^/{NATR) and to nidf^'^d''^^'^/{NATR). We studied the sizes of the Casimir 
scaling violation for various representations of the static sources for G = SU (3). We find 
that they are fairly suppressed and fully consistent with the present bounds from lattice 
calculations, even if we take into account strong scale-dependences of the predictions. 
(Through a diagrammatic analysis, we observe that the suppression of the Casimir scaling 
violation is related to the charge conjugation symmetry.) 

We also discussed possibilities that our predictions for the effects of the Casimir scal- 
ing violation may be tested for representations with large dimensions or in theories with 
larger gauge groups. If these tests are unreachable in near future, one may alternatively 
perform an indirect test of the contribution of the d'^'^'^d^'^'^ / (NaTr) term, using the 
potentials for various representations of large gauge groups. If the perturbative pre- 
diction is excluded by lattice computations in these tests, it means that (at least) for 
the Casimir scaling violating effects non-perturbative contributions are dominant in the 
relevant distance region. Oppositely, suppose the perturbative prediction is positively 
confirmed, for instance, concerning the proportionality of the effects to the specific color 
factors. In this case, the tiny Casimir scaling violation may serve as another discriminant 
of models, which attempt to explain the nature of QCD in the long-distance as well as 
in the intermediate-distance regions. 

It is expected that the violation of the Casimir scaling of the potential has some 
connection with the string breaking phenomenon, since the Casimir scaling property is 
known to be violated in this phenomenon. For this reason, it may be meaningful to 
focus on the Casimir scaling violating effects in the aforementioned class of models, even 
if they are tiny at r < 1 fm. In fact, the string breaking phenomenon is among the 
major dynamical issues of QCD between the intermediate and long-distance regions to 
be explained by these models. We may entertain some speculation from the perturbative 
computation. The diagrams which contribute to the Casimir scaling violations include 
those which (naively) are relevant in the string breaking phenomenon (for instance, 
the first diagram of Fig. [7] in the case R = F, and the second diagram in the case 
R = A). Moreover, the effect of the Casimir scaling violation is predicted to decrease 
the force between the static charges, indicating that screening of these charges is taking 
place. Although it is beyond the scope of perturbative QCD to treat the string breaking 
phenomenon, some hints may be revealed by investigating further the Casimir scaling 
violation in the distance region where perturbative treatment is still valid. 
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